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ABSTRACT 

We present the results of a 430-MHz survey for pulsars conducted during the upgrade 
, to the 305-m Arecibo radio telescope. Our survey covered a total of 1147 deg 2 of sky 

using a drift-scan technique. We detected 33 pulsars, 10 of which were not known prior 
to the survey observations. The highlight of the new discoveries is PSR J0407+1607, 
which has a spin period of 25.7 ms, a characteristic age of 1.5 Gyr and is in a 1.8-yr 
orbit about a low-mass (> 0.2 Mq) companion. The long orbital period and small 
\ eccentricity (e = 0.0009) make the binary system an important new addition to the 

ensemble of binary pulsars suitable to test for violations of the strong equivalence 
principle. We also report on our initially unsuccessful attempts to detect optically the 
J> , companion to J0407+1607 which imply that its absolute visual magnitude is > 12.1. 

If, as expected on evolutionary grounds, the companion is an He white dwarf, our 
' non-detection imples a cooling age of least 1 Gyr. 
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1 INTRODUCTION 

Recycled pulsars — neutron stars that are thought to have 
been spun up by the accretion of matter from a binary com- 
panion — are interesting objects to study. In addition to 
their application as excellent clocks for high-precision tim- 
ing experiments (Davis et al. 1985), progress in our under- 
standing of their origin and evolution continues to be made 
by large-area surveys of the radio sky. 

Following Wolszczan's (1991) discovery of two recy- 
cled pulsars at high Galactic latitudes, Johnston & Bailes 
(1991) demonstrated that the local observable population 
of these objects should be isotropic, and that low-frequency 
(~ 400 MHz) large-area surveys would be an excellent means 
of probing this population. This prediction was largely borne 
out by the discovery of ~ 50 recycled pulsars from a number 
of large-area surveys undertaken in this frequency range by 
many groups around the world (see, e.g., Camilo 1995 for a 
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review). More recently, sensitive surveys at higher radio fre- 
quencies (~ 1.4 GHz) continue to discover millisecond and 
recycled pulsars in large numbers at intermediate and low 
Galactic latitudes (Edwards et al. 2001; Faulkner et al. 2004) 
where interstellar propagation effects limit the Galactic vol- 
ume surveyed at lower frequencies. 

The 305-m Arecibo telescope has played a key role in 
low-frequency surveys. Between 1993 and 1998, the upgrade 
to the telescope meant that it was often parked at a fixed 
azimuth for extended periods. This provided a unique oppor- 
tunity to survey the sky as it drifts through the stationary 
telescope beam at the sidereal rate. At 430 MHz, the sky 
coverage for a declination S is 60 cos 5 deg 2 per day. As a 
number of groups were interested in drift-scan surveys (Fos- 
ter et al. 1995; Camilo et al. 1996; Lommen et al. 2000; 
McLaughlin et al. 2002), the available sky was divided up 
into five zones and the zenith angle of the telescope was 
changed daily to cover each region uniformly. 

In this paper, we report on an analysis of data in 
the regions surveyed by the "Space Telescope Science In- 
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Table 1. Declination strips covered in the STScI/NAIC survey. 
For each strip we list the central zenith angle (ZA), gain (G), 
system temperature (T sys ) of the 430-MHz line feed, approximate 
limiting flux density (5 m i n ), area of sky covered and the number 
of new and previously known pulsars detected (new+old). 



S range 
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(deg) 
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(deg 2 ) 
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17 
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0.6 
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0.5 
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0+5 
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12.5 
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0+1 


36-37 


18 
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1.0 


11 
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stitute and National Astronomy and Ionospheric Center" 
(STScI/NAIC) collaboration. A preliminary account of this 
work was given by Xilouris et al. (2000) . In Section|2] we de- 
scribe the survey observations and data analysis procedure. 
The results, follow-up timing and implications of the ten 
new pulsars found in the survey, including the new binary 
pulsar J0407+1607, are detailed in Sections [3] and [I] 



2 SURVEY OBSERVATIONS AND RESULTS 

The survey observations were carried out between April 
1994 and March 1995 using the 430-MHz carriage house line 
feed. While the nominal sensitivity of this receiver for obser- 
vations close to the zenith is 16 K Jy _1 , these observations 
were made at a variety of zenith angles and during periods 
when the feed was slightly out of focus due to the mainte- 
nance work. As a result, the gain of the telescope varied 
considerably in the range 10-16 K Jy _1 . In drift-scan mode, 
a source passes through the 3 dB points of the primary 
beam (10 arcmin) in about 40 s. The declination ranges 
assigned to the STScI/NAIC collaboration, area coverage 
and approximate system parameters are summarised in 
Table The sky coverage is shown in Fig.0 An ASCII file 
containing the exact sky positions observed is available at 
http:/ /www. blackwellpublishing.com/products/journals/ 
suppmat /MNR/MNR9005 /MNR9005sm.htm. 

The incoming radio signals were sampled by the line 
feed as two orthogonal circular polarisations which, af- 
ter amplification and down-conversion, were passed to a 
32 x 250-kHz analogue filterbank spectrometer which sums 
the two polarisation channels prior to detection with an RC 
time constant of 333 [is. The resulting total-power signals 
were sampled every 250 fis before being digitised with 3-bit 
precision and written to magnetic tape for off-line analysis. 

The data were processed in overlapping 33-s segments 
using a derivative of the software package developed for 
earlier Arecibo pulsar searches (Fruchter 1989; Nice 1992; 
Nice, Fruchter & Taylor 1995). In brief, the data were de- 
dispersed for 448 different assumed values of dispersion mea- 
sure (DM) in the range 0-504 cm pc . The resulting time 
series were then Fourier transformed so that the amplitude 
spectra could be searched for periodic signals buried in the 
noise. To increase sensitivity to signals with narrow duty cy- 
cles (i.e. many harmonics in the Fourier domain) the spec- 



Table 2. STScI/NAIC survey detection statistics. Listed are 
the 10 newly discovered pulsars (marked with asterisks), 23 pre- 
viously known pulsars detected and 6 previously known pulsars 
which lie within 5 arcmin of a pointing which were not detected. 
For each pulsar, we give the period (P), DM, 400-MHz flux den- 
sity (S400) i offset (A) from the beam centre, observed S/N and 
the predicted S/N. Blank S/N entries show the non detections. 
Flux densities marked with a f have been estimated from cata- 
logued S1400 values assuming a spectral index of -1.6. 
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tra were incoherently summed so that up to 2, 3, 4, 8 and 
16 harmonics were combined. Signal-to-noise ratios (S/Ns) 
were computed from the amplitude of each point in the spec- 
trum and the local root-mean-square deviation. The Fourier 
components of the 100 strongest events were then inverse 
transformed back into the time domain to compute the S/N 
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Figure 1. Survey sky coverage equatorial coordinates showing new discoveries (rilled triangles) previously known pulsars detected by 
the survey (circles) and previously known pulsars not detected (squares). Dotted lines show constant Galactic latitudes of ±20°. 



of the resulting pulse profile. The profiles, periods, S/Ns and 
best DMs of each candidate were saved for visual inspection. 

The sensitivity of the survey was dependent on the tele- 
scope's zenith angle and hence varied for the different decli- 
nation strips. In Table Q we quote our theoretical minimum 
detectable flux density (S'min) to long-period pulsars with 
duty cycles of 4% using the radiometer equation (see, e.g., 
Lorimer & Kramer 2005). This calculation assumes a sky 
background temperature of 30 K and that the pulsar tran- 
sits through the beam centre. These values should therefore 
be taken as lower limits to the true survey sensitivity which 
degrades for larger background temperatures, pulse widths 
and shorter pulse periods. Typical values for millisecond pul- 
sars, for example, are expected to be > 3 mjy. A more de- 
tailed discussion on the sensitivity of the observing system 
can be found in Camilo, Nice & Taylor (1996). 

The data analysis resulted in the detection of 33 pulsars 
in 1147 deg 2 of sky. Ten of these pulsars were previously 
unknown prior to the survey observations. In Table [5] we 
list the S/Ns of the pulsars along with an estimate of the 
theoretical S/N based on the nominal survey parameters 
and accounting for position offsets from the centre of the 
(assumed Gaussian) telescope beam power pattern. In gen- 
eral, for S/Ns below 20, the observed and theoretical values 
are in reasonable agreement. For some of the bright pulsars 
(e.g. B0823+06), the detected S/Ns are much less than the 
theoretical values. We are unsure of the reasons for this, and, 
unfortunately, the raw data are no longer available. Given 
our short integration times, it is possible that the observa- 
tions were affected by adverse scintillation and/or nulling. 

Also included in Table [5] are 6 previously known pul- 
sars that lie within 5 arcmin of the searched area but were 



not detected. Although one of these pulsars, J1909+0616, 
made the list of candidates with a S/N = 7.6, this was below 
the nominal survey threshold S/N of 9 and does therefore 
not count as an blind detection. Only two of the remaining 
five pulsars have theoretical S/Ns above this threshold. For 
PSR J1920+1110, we are unsure of its 430-MHz flux den- 
sity which has been extrapolated from the 1400-MHz value 
(Morris et al. 2002) assuming a spectral index of -1.6. As this 
pulsar was not detected in more sensitive 430-MHz obser- 
vations by Bhat et al. (2004), either the spectrum is flatter 
than assumed or the pulse is significantly affected by scat- 
tering at 430 MHz. The remaining pulsar which is nominally 
above our detection threshold, J1951+1123, has a period of 
5.1 s and is known to null for several periods (Nice 1992). 
This is the most likely explanation for our non-detection. 

In summary, given the uncertainties in pointing and un- 
known variations in received focus during these observations, 
this simple analysis of the detection statistics shows that our 
survey sensitivity is consistent with that reported by other 
groups (see, e.g., Camilo et al. 1996). 



3 FOLLOW-UP OBSERVATIONS 

Confirmation and follow-up observations of the new discov- 
eries began at end of the Arecibo upgrade in 1998. These 
made use of the Penn State Pulsar Machine (PSPM; Cadwell 
1997), a 128x60-kHz channel analogue filterbank spectrome- 
ter which records either continuously-sampled data ( "search 
mode") or synchronously-averaged pulse profiles at the pre- 
dicted pulse period ( "timing mode" ) . 

As in the original survey observations, we used the 430- 
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Table 3. Parameters for the nine isolated pulsars derived from the TEMPO timing analysis. Figures in parentheses are l-cr uncertainties 
in the least-significant digits as reported by TEMPO. The arrival times from which these ephemerides are derived are freely available 
on-line as part of the European Pulsar Network (EPN) database (http://www.jb.man.ac.uk/~pulsar/Resources/epnl. 
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Figure 2. Integrated 430-MHz pulse profiles for nine of the newly 
discovered pulsars derived from a phase-coherent addition of the 
individual PSPM observations. Each profile shows 360 degrees of 
rotational phase. These profiles are freely available as part of the 
EPN database. A pulse profile for PSR J0609+2130 has already 
been published by Lorimer et al. (2004) and is not included here. 



MHz line-feed receiver for the bulk of our follow-up work. 
Following initial confirmation observations which verified 
the existence of each pulsar in search mode, we generated 
a preliminary ephemeris which was subsequently used to 
obtain timing-mode data. The basic observing and anal- 
ysis procedures used were identical to those described by 
Lorimer, Camilo & Xilouris (2002). 

In order to determine accurate spin and astrometric pa- 
rameters for each pulsar, we carried out a dedicated cam- 
paign of observations using the PSPM in timing mode. For 
each observation the set of 128 folded pulse profiles (one 
profile per frequency channel) was appropriately delayed 
with respect to the centre observing frequency to account 
for the effects of interstellar dispersion. Standard software 
tools (Lorimer 2001) implemented this procedure to form 
a de-dispersed time-tagged integrated pulse profile for each 
observation. A pulse time of arrival (TOA) for each pro- 
file was then obtained by cross correlating it in the Fourier 
domain with a high S/N template profile formed from the 



sum of many individual observations (Taylor 1992). The fi- 
nal template profiles are shown in Fig. 

The TO As obtained for each pulsar were transformed to 
equivalent arrival times at the solar system barycentre and 
fit to a simple spin-down model using the TEMPO 1 soft- 
ware package. After several iterations, an initial ephemeris 
was obtained which was then used to form an improved tem- 
plate profile from the sum of each pulse profile appropriately 
phase-aligned according to the timing model. The resulting 
high S/N template was then used to obtain TO As of higher 
precision which were analysed with TEMPO in a further it- 
eration to produce the final timing solution. Refined values 
of DM for each pulsar were obtained by forming four 1.92- 
MHz-wide sub-bands from each PSPM profile. The DM was 
derived from a fit to the TO As from each sub-band keeping 
the spin and astrometric parameters fixed at the final tim- 
ing solution values. For nine of the new pulsars, the TOAs 
were adequately fit using an isolated pulsar spin-down model 
and the resulting parameters from each fit are presented in 
Table El 

In the time since the original survey observations 
and follow-up, five of the new pulsars have been inde- 
pendently discovered in other surveys. The binary pulsar 
J0407+1607 discussed below was detected during a survey 
of gamma-ray error boxes at Parkes (Roberts et al. 2004). 
PSRs J0609+2130 (Lorimer et al. 2004) and J1549+2113 
(Foster et al. 1995; Lewandowski et al. 2004) were discovered 
in contemporary 430-MHz Arecibo drift-scan surveys. Fi- 
nally, PSRs J1848+0604 (Faulkner, private communication) 
J1905+0616 (Morris et al. 2002) were detected by the Parkes 
multibeam survey. PSRs J1549+2113 and J1848+0604 have 
no previously published timing solution and are presented 
for the first time in Table El For PSR J0609+2130, we have 
added an additional year of PSPM TOAs to those already 
published by Lorimer et al. (2004) . The longer time baseline 
of the resulting ephemeris presented in Table El has signif- 
icantly improved the precision of the derived parameters. 
For PSR J1905+0616, we have included Jodrell Bank TOAs 
collected at 1412 MHz by Morris et al. (2002) in our final 
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Table 4. Measured parameters for PSR J0407+1607 from 71 
TOAs from MJDs 52304 (January 30, 2002) to 53292 (October 
14, 2004). The numbers in parentheses are the 1-a TEMPO un- 
certainties in the least significant digit quoted. 



Parameter 


Value 


Right ascension (J2000) (h:m:s) 


04:07:54.939(3) 


Declination (J2000) (°:':") 


16:07:16.4(2) 


Spin period, P (ms) 


25.70173919463(2) 


Period derivative, P (xlO -20 ) 


7.9(3) 


Epoch of period (MJD) 


52799 


Dispersion measure, DM (cm~ 3 pc) 


35.65(2) 


Projected semi-major axis, x (s) 


106.45026(2) 


Orbital period, Pi, (days) 


669.0704(1) 


Orbital eccentricity, e 


0.0009368(6) 


Longitude of periastron, u> (deg) 


291.74(2) 


Epoch of periastron, Tq (MJD) 


52774.69(3) 



timing analysis. This has resulted in an improvement on the 
ephemeris published by Morris et al. 

The most interesting of the new discoveries is 
PSR J0407+1607. The intermediate 25.7-ms spin period of 
this pulsar suggested that it was either a young or recycled 
object. Following initially unsuccessful attempts to fit an iso- 
lated spin-down model to the data, sets of closely-separated 
TOAs were used to obtain estimates of the barycentric pulse 
period for a given epoch. These periods were then fit to 
a simple Keplerian model which revealed 1.8-yr sinusoidal 
variations implying that the pulsar was a member of a long- 
period circular orbit binary system. 

While our early PSPM timing-mode observations could 
be used to obtain a preliminary orbital ephemeris, they were 
not of sufficient quality or density to obtain a phase-coherent 
timing solution. In January 2002, we began a regular observ- 
ing campaign to collect data with the PSPM in search mode. 
These data were then folded off-line to produce a set of pulse 
profiles and hence TOAs for a timing analysis. The advan- 
tage of this approach is that the data could be re-folded later 
when the precision of the orbital parameters improved. 

Despite this improved approach, small uncertainties in 
the assumed position resulted in significant covariances be- 
tween the orbital parameters and position fits in TEMPO. 
Fortunately, as mentioned above, this pulsar was indepen- 
dently discovered by Roberts et al. (2004). Follow-up 1400- 
MHz Arecibo observations from that survey resulted in a 
positional determination which was used as a starting point 
for the final timing analysis. The final phase-coherent tim- 
ing solution presented in Table [3] is the result of a simulta- 
neous fit for spin, astrometric and binary parameters. The 
observed minus computed TOA residuals for this solution 
are featureless with a root-mean-square value of 16 fis. The 
DM was determined by including 1400-MHz observations 
of the pulsar with the Wideband Arecibo Pulsar Processor 
(Dowd, Sisk & Hagen 2000) and, as for the other pulsars, 
keeping all other model parameters constant in the TEMPO 
fit. 

To obtain an estimate of the transverse speed of 
J0407+1607, we performed a scintillation analysis on an ob- 
servation at 430 MHz using the Arecibo Observatory Fourier 
Transform Machine (AOFTM), a fast-dump digital spec- 
trometer which provides 1024 frequency channels across a 




Pulse Phase 



Figure 3. Upper panel: full Stokes profile for PSR J0407+1607 
at 430 MHz showing total intensity (solid line), linear polarisation 
(dashed) and circular polarisation (dotted) as a function of pulse 
phase. Lower panel: position angle versus pulse phase. 



10-MHz bandwidth (see, e.g., Mclaughlin et al. 2002). Per- 
forming a standard autocorrelation analysis on the pulsed 
intensity as a function of time and radio frequency (for de- 
tails of this technique see Lorimer & Kramer 2005) we find 
the characteristic timescale and bandwidth of scintillation 
to be 843 ± 6 s and 190 ± 2 kHz respectively. Adopting the 
expressions given by Cordes & Rickett (1998) for a uniform 
scattering medium, and taking a distance of 1.3 ± 0.2 kpc 
derived from the DM and the Cordes & Lazio (2002) elec- 
tron density model, we estimate the transverse speed of 
PSR J0407+1607 to be 35 ± 5 km s _1 . Although this re- 
sult is model dependent, it implies a transverse speed which 
lies at the low end of the distribution known for millisecond 
pulsars (Nice & Taylor 1995). In future it should be possi- 
ble to measure the proper motion of this pulsar directly by 
extending the existing baseline of timing observations. The 
predicted proper motion at this distance is 6 ± 1 mas yr _1 . 
Through simulations with TEMPO, we estimate that a di- 
rect proper motion measurement of this magnitude through 
timing techniques will be possible within 3 yr. The effect of 
proper motion on the observed semi-major axis should also 
be measurable on a similar timescale, allowing us to place 
constraints on the orbital inclination angle in future. 

The PSPM and AOFTM data do not provide 
any information about the polarisation properties of 
PSR J0407+1607. To investigate this, in March 2002, we 
observed the pulsar for 78 min at 430 MHz using the Prince- 
ton Mark IV data recording system (Stairs et al. 2000) in 
narrow-band (5 MHz) mode. The Mark IV data were co- 
herently de-dispersed off-line and the data were calibrated 
using observations with a noise diode of known flux den- 
sity and polarisation. The total intensity, linear and circu- 
larly polarised pulse profiles are shown in Fig. [3] Relative to 
the total intensity, we find that both the linear polarisation 
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Table 5. Observed and derived parameters for the ten pulsars discovered in the survey. For each pulsar, we list the Galactic longitude 
(/) and latitude (fc), the flux density at 430 MHz (S430), the pulse width measured at 50 and 10% of the peak intensity (respectively W50 
and W10), the equivalent width of a top-hat pulse having the same area and peak intensity of the profile (W eq ), the distance derived 
from the dispersion measure assuming the Cordes & Lazio (2002) electron density model (D), as well as the base-10 logarithms of the 
430-MHz luminosity (L 430 = 5430-D 2 ), the characteristic age (r = P/(2P)), the surface dipole magnetic field (B = 3.2 X 10 19 (PP) 1 / 2 G) 
and the rate of loss of spin-down energy (E = 3.95 X 10 46 P/P 3 erg s _1 ). 
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and modulus of the circularly polarised component have an 
average value of 13%, typical of other millisecond pulsars 
(Xilouris et al. 1998; Stairs, Thorsett & Camilo 1999). Also 
shown in Fig.^is the position angle of the linear polarisation 
as a function of pulse phase. As for many other millisecond 
pulsars the position angle does not follow a smooth S-shape 
variation expected from the rotating-vector model of Rad- 
hakrishnan & Cooke (1969). A determination of the emission 
geometry of this pulsar is, therefore, non-trivial. 



4 DISCUSSION 

Various observed and derived parameters of the new pulsars 
are summarised in Table [3] where we also include estimates 
of each pulsar's 430-MHz flux density. These were derived 
from the PSPM profiles using identical procedures to those 
described by Lorimer et al. (2002). 

4.1 General remarks 

In common with results from other Arecibo drift-scan sur- 
veys (Foster et al. 1995; Camilo et al. 1996; Lommen et 
al. 2000; McLaughlin et al. 2002), the new pulsars are rela- 
tively faint, nearby objects. Their inclusion in future popula- 
tion syntheses should provide valuable constraints on the low 
end of the radio pulsar luminosity function. In 1147 deg 2 of 
sky searched, we detected five recycled pulsars: J0407+1607, 
J0609+2130, B1534+12, B1913+16 and B1937+21. Of the 
two new discoveries, the binary pulsar J0407+1607 is dis- 
cussed below, while the isolated pulsar J0609+2130 was al- 
ready discussed in detail by Lorimer et al. (2004) . We simply 
note here that our detection rate of one recycled pulsar per 
~ 230 deg 2 is consistent with that of other Arecibo surveys 
(see Camilo 1995). 

The youngest pulsar in our sample, with a characteristic 
age of 116 kyr is PSR J1905+0616. As noted and discussed 
by Kramer et al. (2003), this object lies within the bounds 
of the EGRET source 3EG J1903+0550. An association is 
unlikely simply due to the large density of pulsars in this 



region and hence high probability of a chance aligmnent. In 
addition, as noted by Kramer et al., the spin-down energy 
loss of J1905+0616 is not large enough to be a likely source 
of the gamma-ray emission seen. 

4.2 PSR J0407+1607 and 3EG J0407+1610 

PSR J0407+1607 is coincident with the unidentified 
EGRET source 3EG J0407+1610 (Hartman et al. 1999), 
a relatively weak gamma-ray source with mean flux F n = 
2.4 x 10 -10 erg cm -2 s _1 above 100 MeV (assuming a pho- 
ton spectral index of -2). While the variability index (see 
McLaughlin et al. 1996 for a definition) is only 0.19, and con- 
sistent with other EGRET sources known to be associated 
with pulsars, our P measurement of J0407+1607 implies a 
rate of loss of spin-down energy E — 1.9 x 10 32 erg s -1 . 
This can be used to test the validity of an association 
with 3EG J0407+1610 by calculating the required efficiency 
rj of conversion of the pulsar's spin-down luminosity into 
gamma rays. At the DM-derived distance of 1.3 kpc, and 
assuming a gamma-ray beaming fraction of 1/47T sr, we find 
an incompatible value of rj = (F y D 2 /E) ~ 20. A much 
smaller distance and/or beaming fraction would be required 
so that 77 < 1. We therefore conclude that PSR J0407+1607 
is most likely not the source of the gamma rays from 
3EG J0407+1610. 

4.3 PSR J0407+1607 and SEP violations 

The orbital parameters of the J0407+1607 binary system 
are consistent with standard models of binary evolution (see, 
e.g., Phinney & Kulkarni 1994) in which the radio pulsar has 
been spun up by the accretion of matter while its companion 
star was in a red giant phase. Indeed, the observed orbital 
eccentricity is in excellent agreement with the prediction 
from Phinney's (1992) relationship between eccentricity e 
and orbital period Pj,. 

Millisecond pulsar-white dwarf binaries like 
PSR J0407+1607 make excellent probes of the strong 
equivalence principle (SEP). The SEP states that free fall 
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of a body is completely independent of its gravitational self 
energy. If the SEP is violated, then the ratio of inertial 
mass to gravitational mass of a falling body differs from 
unity by an amount A. As originally proposed by Damour 
& Schafer (1991), and by analogue with Nordvedt's earlier 
tests of the SEP in the Earth-Moon system (Nordvedt 
1968a,b), a violation of the SEP for a pulsar-white dwarf 
binary would mean that both objects "fall" in the local 
Galactic gravitational field at a different rate. As a result, 
a violation would cause a polarisation of the orbital plane 
toward and parallel to the centre of the Galaxy. 

Following Damour & Schafer's early work, Wex (1997, 
2000) derived limits on A from an ensemble of binary pulsar- 
white dwarf systems. As the figure-of-merit for this cal- 
culation is Pb/e 2 , long-period circular orbit binaries like 
J0407+1607 are particularly important members of the en- 
semble used in the test. Recently, a robust Bayesian analysis 
by Stairs et al. (2005) which uses all known binary pulsars 
relevant for this test including J0407+1607 find a limit on 
SEP violations of A < 5.5 x 10" 3 . This new limit is com- 
parable to the weak-field limits found in the Earth-Moon 
system (Nordvedt 1968a, b). 



4.4 The companion to PSR J0407+1607 

The observed mass function of PSR J0407+1607 is 2.89 x 
10 -3 M©. Assuming a mass for the pulsar in the range 1.3- 
1.6 Mq , the implied companion mass is 0.2 Mo /sin i, where i 
is the unknown inclination angle between the plane of the or- 
bit and the plane of the sky. From extensive studies of other 
millisecond pulsar binaries with similar orbital characteris- 
tics to J0407+1607 (for a recent review, see van Kerkwijk et 
al. 2005) it is most likely that the companion to J0407+1607 
is a low-mass white dwarf. An inspection of archival Palomar 
data in the digitized sky survey archive 2 reveals no optical 
counterpart down to the plate limit (approximately 19th 
magnitude), with the closest object in the field being 8 arc- 
sec (50cr) offset from the radio position given in Table 3] 

To search for a fainter optical counterpart, the field sur- 
rounding PSR J0407+1607 was observed on October 2003 
with the 2.3-m Bok telescope at Kitt Peak using the Stew- 
ard Observatory imaging camera in combination with the 
Harris B and V filters. Four 27-min B and four 20-min V 
images were obtained under non-ideal seeing (~ 1.5 arcsec) 
conditions. All images were bias-subtracted and flat-fielded 
using the IRAF package 3 . For photometry calibration, we 
selected all stars in a 3 x 3 arcmin box around the pulsar 
location from the USNO-B1.0 Catalog 4 . 

The detection limit was determined to be magnitude 21 
in B. This is in good agreement with the 20.9 magnitude of 
the star identified as 1061-0042921 in the USNO-B1.0 cata- 
logue detected at the noise level in our images. Four 20-min 
exposures were aligned and added to produce the final V im- 
age. We used the published visual magnitude of the nearby 
stars SA 94-171 and GD 71 (Landolt 1992) to determine 
the relationship between instrumental and aperture magni- 



2 http://www-gsss.stsci.edu/DSS 

3 http://iraf.noao.edu 

4 http:/ /www. nofs.navy.mil/data/fchpix 



tudes. From this procedure, we estimate a limiting apparent 
visual magnitude for detection of 23. 

At the nominal distance of 1.3 ± 0.2 kpc, our upper 
limit corresponds to an absolute visual magnitude of at least 
12.1. Assuming the companion to be an He white dwarf, we 
can use the cooling curves computed by Hansen & Phinney 
(1998) to place a lower limit on the age of the white dwarf 
companion. From an inspection of Fig. 13 in their paper, 
our current non-detection constrains the cooling age to be 
;S 1 Gyr, compatible with the 1.5 Gyr characteristic age of 
J0407+1607. More stringent constraints should be available 
soon following deeper observations of this system now un- 
derway with the Palomar telescope which will be described 
elsewhere. 
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